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SUMMARY 


BRL has undertaken a program to investigate the adequacy of methods 
currently available for predicting the effects of topographic features on 
the blast wave from nuclear weapons. 
As part of that program, under Contract 
No. PMPDG-77-C-0721 , Scientific Service, Inc. (SS!) was requested to aid in 
the design of a model based on real terrain, and to make predictions of the 
effects of this model terrain on blast waves generated by 454 gm high ex- 
plosive (HE) charges. 
The size of the HE charges, and their height of burst 
over the model lO.6 m} were chosen to simulate the burst of a 125 kt weapon 
at a height of 300 m over prototype terrain. 


The work undertaken by SSI included an extensive review of the more 
than 15 years of research work on which the current prediction methods are 
based. 
Thi s revi ew, and the s'lIbsequent appl i cati on of currently avail able 
prediction techniques, led to the following conclusions: 


• There are unresolved differences between two techniques for predict- 
ing shock front overpressures as affected by terrain, both of which at one 
time or another had been recommended for field use. 


• Certain terrain features can generate transients of very short dura- 
tion which place severe requirements on the time response of the systems for 
measuring the blast pulses from small HE charges. 
(These transients can be 
important in full-scale situations.) 


• 
No techniques were available for predicting the effects of terrain 
on pressure behind the shock front in complex terrain. A new method was 
devised to do this. 


• 
In certain situations. maximum pressures can occur behind the front 
of the blast wave, and these can differ substantially from shock front pres- 
sures, which are usually thought to be maximum. 


• Errors in the current prediction curves were found and corrected. 
Because of the severe requirements placed on the pressure recording systems, 
the model study may not answer all unresolved questions. 
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I. 
INTRODUCTION 


The characteristics of a blast wave propagating along the ground surface 
change when it encounters a topographic feature: 
pressures in the pulse gen- 
erally increase when the wave encounters a rising slope or proceeds down the 
axis of a valley with parallel or converging sides; pressures generally decreasE 
when the wave encounters a falling slope. 
Under certain conditions, these 
changes in pressure can be significant. 


Techniques were developed some years ago for predicting these terrain- 
related changes, and a few tests confirming their accuracy were conducted. 
Recently, however, questions arose about the adequacy of the techniques, that 
is, whether they could deal with the practical problems that might be encounter- 
ed in the field. 


For this reason, BRL has undertaken a program aimed at evaluating the 
adequacy of currently available techniques. 
The basic plan for this evalua- 
tion is quite straightforward. Predictions would first be made of the character 
istics of blast waves from 454 gm charges exploded over a small scale model con- 
taining topographic features resembling those of real terrain. 
(Burst location~ 


and heights would be based on a prototype full-scale situation involving weapon 
yields of 125 kt and burst heights of 300 m.) 
Tests would then be conducted, 
and the experimental results would be compared with the predictions. 
Both pre- 
di cti ons and experimental measurements woul d be restri cted to bl ast overpres- 
sures. 


As part of this program, BRL awarded Contract DAAD05-77-C-0721 toScientifi 
Service, Inc. (SSI) to: 
design the model; select locations for the pressure 
transducers in the model; make recommendations about instrumentation; and, 
finally, to make predictions of complete pressure·time histories for the burst 
conditions and transducer locations selected. This report summarizes these 
activities. It describes the techniques now available to make predictions of 
blast wave pressure changes caused by terrain features; discusses the reasons 
for choosing a particular model and measuring station layout; and applies the 
prediction techniques to the particular situtations to be tested. 


II. 
PREDICTION METHODS 


A. 
Background 


Methods currently available for predicting the effects of terrain on air 
blast stem from programs, begun in the early 1950 
l s and extending for some 15 
years, involving theoretical analyses, and experiments with nuclear weapons, 
shock tubes and small high explosive charges. 


The program was initiated after some terrain effects were noted during 
early weapons tests (Ref. 1) and experiments using simple (constant slope) 
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terrain models and small high explosive charges, Vfere conducted to investigate 
the magnitude of the effects (Refs. 2, 3,4). The information from these ex- 
periments was used to establish methods for predicting changes in overpressure 
at and behind the front of shock waves as a function of incident overpressure, 
slope angle (the angle of steepest ascent ), and azimuth angle (the angle between 
the direction of shock propagation and a line drawn perpendicular to the slope 
contours) . 


These methods were then used to make predictions of effect in two complex 
terrain models (Ref. 5) and in real terrain (Ref. 6). There was close corre- 
lation between measured and predicted results in the terrain models. 
On the 
full scale test on real terrain) rising slope terrain effects were masked by 
the presence of precurser phenomena, in whi ch i nteracti on between the thermal 
pulse from the burst and the ground surface caused the blast wave to develop 
an irregular configuration generally characterized by a rounded front and by 
higher than normal dynamic pressures. 
On those terrain features (mostly fall- 
ing slopes) which did not see the thermal pulse, predictions compared well 
with measurements (Ref. 7). 


Subsequent to the full-scale test, theoretical work was undertaken to 
develop techniques for predicting terrain effects on dynamic pressure (Refs. 8, 
9). 
As a by-product of this work, an analytic approach to predicting shock 
front overpressure (Ref. 10) was found to agree remarkably well with the results 
of an earlier extensive series of shock tube measurements (Ref. 11). 
Because 
of its sound theoretical basis and good agreement with measurement, a predic- 
tion technique based on this approach supplanted the earlier one based on 
meas urements from tests wi th small charges. 
Generally the two techniques 
yield essentially the same results at higher incident overpressures, but at 
lower incident overpressures the analytic technique predicts higher pressures 
on rising slopes, and lower pressures on falling slopes than does the technique 
based on small-charge tests. 


The special problem of channeling -- in which a shock wave propagates 
along the axis of a valley with parallel or converging sides -- was then studied 
both theoretically and experimentally with charges weighing 110 gm, and a method 
for predicting channeling effects was developed (Ref 12). 
Finally, methods for 
dealing with a variety of real terrain features were developed (Ref. 13), and 
now form the basis for the techniques suggested for general use (Ref 14). 


B. 
General Considerations 


When a blast wave propagating over the ground surface encounters a negative 
or falling slope, it undergoes a phenomenon known as diffraction.* The wave 


* In general discussions of gas dynamics, the term "diffraction" is most fre- 
quently used to refer to effects of any slope change, negative or positive. 
Less confusion results if the term is used to describe effects of negative 
slope changes only, and the term "reflection" used to describe effects of 
positive slope changes. 
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expands to stay in contact with the slope, becomes curved, and pressures at and 
behind the shock front decrease. A typical shock pattern after diffraction is 
sketched in Fig. lA. 
If the blast wave encounters a positive or rising slope, 
it undergoes reflection*, which can be of two types. 
A steep slope causes 
"regular reflection" (also known as "two-shock" reflection), in which a single 
reflected shock forms behind the incident shock as shm·Jn in Fig. lB. 
A shallow- 
er slope causes "Mach reflection ll in which a third shock, the ~ilach stem, forms 
between the ground surface and the conjunction of incident and reflected shocks 
producing a pattern like that of Fig. Ie. 
The demarcation between conditions 
that cause regular Mach reflection is not sharp; indeed, there are certain con- 
ditions under which both types of reflection can occur. 
In either case, howevE 
pressures in the shock wave at the ground surfaces are greater than they are in 
the incident shock wave. 


Hhen a blast wave propagates down the axis of a valley \'iith either parallE 
or converging sides, (that is, undergoes channeling) it is confined by the 
valley walls and pressures within the wave become greater than they would be 
if it were propagating over flat terrain. 


The changes a wave undergoes when it encounters a slope change are depende 
on the strength of the shock (i.e., overpressure behind the shock front rela- 
tive to ambient pressure ahead of the front), and the angle of the slope change 
or slope angle.* 
This angle would be the angle es in Fig. 1, if the direction 
of propagation of the shock in the three sketches of Fig. 1 were normal to the 
slope contours. 
In the general case, where the direction of propagation is not 
normal to the slope contours, an "effective slope angle" e is used, defined as 


sin e = (sin 8s) (cos cp) 
(la 


where cp is the angle between the direction of shock propagation and the normal 
to the slope contours. 
The angle cp is shown in the inset to Fig. 2, a plot 
relating e, es ' and cp. 


In the follm'iing, prediction graphs for various types of interaction are 
given. 
The most common shock-terrain interactions are those involving Mach re- 
flection and diffraction. These are discussed first, and are followed by more 
limited discussions of regular reflection and channeling. 
The section closes 
with material describing how the prediction graphs are to be used, that is, how 
to make predictions of topographic effects. 


C. 
Prediction Graphs 


1. 
Mach Reflection and Diffraction. 
Methods in current use for predictin 


*' 
Alternatively, the angle between a normal to the shock front and a normal t 
the slope contours, termed the lIangle of incidence ll is frequently used when 
discussing shock waves encountering slopes. This angle equals 90° minus th 
slope angle. 
For discussing topographic effects, the slope angle is more c 
venient to use since it can be derived directly from topographic maps. 
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Fig. 2. 
Effective Slope Angle 0 as a Function of Slope Angle Os for Various Azimuth Angles ~. 


the effects of Mach reflection and diffraction are based on an approximate theory 
developed by Whitham (Ref. 10). 
In that theory it is assumed that the paths of 
air particles remain perpendicular to the shock wave even when it is being dis- 
torted by reflection or refraction. 
t~hile the Whitham theory is thus only approx- 
imate, the excellent agreement between its predictions of shock strengths and 
White's measurements of shock strength in a shock tube (Ref. 11), shown in Fig. 3, 
indicate that the theory has wide applicability. 


The theory involves a function w defined by 


M 
:f - 
1/2 


W = {(M2-\)K} 
dt~ 


1 


where 
M = the t~ach number of the shock 
= the shock propagation velocity/speed of sound 


K = 2{[1 + (y ~ 1) (1 ~ jJ2)] [ 211 + 1 + M- 
2]} -1 


= (y- 1) t12 + 2 _ 
II 


2y t~2 - 
(y - 1) 


(2) 


and 
y = the ratio of specific heats of a gas at constant pressure and volume. 
For air, y = 1.4 


The Whitham function w was integrated by numerical methods especially for 
this program and is shown in Fig. 4, along with the Mach number M, both as 
functions of shock wave pressure.* The change in shock front pressure caused by 
a change in slope angle 8 is derived by finding the value of w for a particular 
value of incident shock pressure Pi from the graph of Fig. 4; adding 8 (in radians) 
to w; then using Fig. 4 again to find a value of pressure corresponding to the 
value of w + 8. 
Note that 8 may be either positive or negative depending on 
whether the shock wave encounters a rising or a falling slope. 


Curves more convenient than those of Fig. 4 for anyone making predictions of 
the effects of topography, are shown in Figs. 5and 6. 
These curves directly relate 
overpressure at the shock front (also called "zero-time" overpressure) on a terrain 


* In this report, the term "pressure" when referring to a shock wave will be 
used interchangeably vdth the term "overpressure", that is, the pressure 
behind a shock front in excess of ambient pressure ahead of the front. 
Thus, 
the "pressure 'l of a weak shock wave approaches zero. 
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5. 
Overpressure at the Shock Front on a Rising Slope (Reflected Overpres- 
sure) vs Peak Incident Overpressure for Shock Waves Undergoing Mach 
Reflection. 
From Whitham Theory. 


Note: 
Numbers in parenthesis identify extreme (limit) slope angles 
along the "Limit for ~lach Ref1ection" line. Other numbers 
on the figure identify the slope angles to which each curve 
applies. 
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Fig. 
6. 
Overpressure at the Shock Front on a Falling Slope (Diffracted Over- 
pressure) vs Peak Incident Overpressure, 
From Whitham Theory. 


Note: 
Numbers on the figure identify the slope angles to which each 
curve applies. 
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feature to peak incident overpressure for various slope angles.* 


Note on Fig. 5, the line labelled "Limit for t~ach Reflection". Generally, 
for a given incident overpressure if the slope angle is greater than the slope 
angle values on the line, regular (two-shock) reflection would occur. 


The prediction curves based on experiments with small HE charges which the 
Whitham curves supplanted are shown in Figs. 7 and 8. 
The Whitham theory curves 
and the small charge curves are in generally good agreement at incident pressure 
levels greater than 70 kPa with rising slopes, and about 150 kPa with falling 
slopes. 
With 
lo\'~er incident pressure levels, the Whitham curves yield higher 
values of pressure on rising slopes, and lower values on falling slopes. 


The small charge studies showed that peak overpressure did not always occur 
at the shock front on a rising slope. This had also been observed in White's 
earlier shock tube work and led to the generation of curves for predicting changes 
in pressure behind the shock front. 
Prediction curves developed for both rising 
and falling slope angles of 10, 20, and 30 deg. are shown in Figs. 9 and 10.** 
Note that the ordinate for these curves is "pressure ratio" (pressure on a slope 
divided by incident pressure) not reflected pressure itself as in Figs. 5 through 
8. 


The abscissas of these curves are values of the ratio t/L, where t is time 
after passage of the shock front over a point of measurement, and L is the dis- 
tance of that point from the location upstream at which the slope changed. 
The 
units (ms/s or s/km) indicate that a shock wave on a slope follows the "self 
similar principle", that is, it grows similar to itself in time from the point 
of slope change. 
As an illustration of the principle, the pressure ratio in a 
shock wave on a slope at a time 0.1 msec after its passage over a point 1 m 
from the start of the slope (t/L = 0.1 ms/m) will be the same as the pressure 
ratio 0.5 msec after passage of the same shock wave over a point 5 m from the 
start of the slope (t/L = 0.5/5 = 0.1 ms/m). 
The principle, which holds strictly 
for two dimensional shock \IJaves with constant pressure behind the shock front, 
was also shown to hold approximately for three dimensional blast waves with de- 
creasing pressures behind the front (Ref. 3), 


Figs. 9 and 10 show directly what the pressure pattern of a two dimensional, 
constant pressure (flat-topped) shock wave vlould look like as a function of time 


* These are similar to the curves given in Ref. 14 with some corrections and 
additions. 
The rising slope (reflected pressure) curves, and the "Limit for 
Mach Reflection 
ll line, required correction at incident overpressure values 
of less than about 20 kPa; the falling slope (diffracted pressure) curves 
required correction over much of their range. 
In addition, for purposes of 
this program, the falling slope curves have been extended to diffracted 
pressure values less than 7 kPa. 


** The curves of Figs. 9 and 10 were derived for equal increments of pressure in 
psi rather than kPa. 
Limits on both time and funds on this project did not 
permit rederivation of the curves for equal pressure increments in kPa. 
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7. 
Overpressure at the Shock Front on a Rising Slope (Reflected Overpres- 
sure) vs Peak Incident Overpressure for Shock Waves Undergoing Mach 
Reflection. 
From Small Charge Studies. 


Note: 
Numbers in parenthesis identify extreme (limit) slope angles 
along the IlLimit for r~ach Reflection 
ll line. Other numbers on 
the figure identify the slope angles to which each curve applies. 
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after it encounters a slope. 
For example, on a 30 deg. rlslng slope (Fig. 9) 
a 14 kPa (2 psi) incident shock wave would become strongly peaked. 
Pressure 


at the front would be about 36 kPa or 5.2 psi (pressure ratio = 2.6) and it 
would decrease by about 50% at tiL of about 0.2 ms/m. 
In contrast, a 690 kPa 
(100 psi) incident wave would have a pressure ratio at the shock front of about 
1. 7. 
This waul d fi rst increase by over 25% at tiL ~ 0.03 ms/m and then stay 
relatively constant. 


Though not apparent, the curves imply that pressure changes caused by a 
sudden slope change spread out in time as the shock wave propagates up a slope. 
(The further up the slope a measurement is made, the later in the pulse will 
a particular effect be noticed.) This phenomenon can substantially alter a 
pressure pulse as it moves over a terrain feature, as can be seen in the three 
examples of Fig. 11. 
The first two examples are for 30 deg. rising slopes; the 
third is for a 30 deg. fall ing slope. 
The basic time scales on the figure are 
for a nuclear burst with a yield of 1 kt, although additional labelling has been 
added to show anticipated effects for nuclear bursts with a yield of 1 Mt, as 
well as high-explosive (TNT) bursts from charges weighing 1 lb (454gm). 


In each example, four pulses are shown. 
The dotted line is the incident 
pulse, i.e., the pulse that would be at the measuring station in the absence of 
a terrain feature. 
For illustrative purposes, peak incident pressure for Fig. 11A 
was chosen as 690 kPa 
(~100 psi); for Figs. lIB and 11C, as 35 kPa 
(~5 psi). 
For 
clarity, (because the high-pressure pulse is of relatively short duration) the 
time scale of Fig. 11A is five times those of Figs. lIB and 11C. 


The incident pulse for Fig. 11A is from Ref. 15; for Figs. lIB and 11C the 
incident pulse shown is calculated from the modified exponential relationship 
given later in this report in the section on predictions. 


The other three curves in each example show the pulses from a burst with a 
yield of 1 kt that would be measured at 1 m, 10 m, and 100 m from the start of 
each slope. 
(For a 1 Mt yield, the distances would be 10 m, 100 m, and 1 km; 
for a 1 lb HE charge they would be 1 cm, 10 cm, and 1 m.) 


In the first case, (rising slope, incident pressure = 690 kPa) the pressure 
peak at the station closest to the start of the slope is behind the shock front 
and is about 20% greater than that at the other two stations. 


In the second case (risingslope, incident pressure = 35 kPa) the values of 
pressure at the peaks (the shock fronts) are the same at the three stations, 
but the pulse at the closest station has a very sharp spike. 


In the last case (falling slope, incident pressure = 35 kPa) three different 
pressure maxima would be observed. 
The maximum at the middle station would be 
about 30% greater than that at the furthest stati on; the maximum at the closest 
station, about 75% greater. 


2. 
Regular Reflection. 
As already noted, for a given incident pressure, if 
the slope angle is greater than the values along the lines labelled "Limit for 
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~1ach Reflectionll in Figs. 5 and 7, regular (two-shock) reflection generally takes 
place. 
(See Fig. 1.) A plot of shock pressure after regular reflection vs inci- 
dent pressure for various values of slope angle is shown in Fig. 12. 
The curves 
were derived from Ref. 16. 
They correct and extend the range of the regular re- 
flection curves of Ref. 14, and are somewhat more convenient to use since reflect- 
ed pressure can be read directly. 


In general, reflected pressures decrease as slope angles increase up to about 
55 deg. 
At incident pressures below about 100 kPa reflected pressures are essen- 
tially constant for all values of slope angles between 55 and 90 deg. 
At inci- 
dent pressures greater than about 100 kPa, however, reflected pressures increase 
as slope angles increase above 55 deg. 
Thus, at incident pressures below about 


100 kPa, the curves labelled 55 deg. and 90 deg. are merged; at incident pressures 
greater than about 100 kPa, the 90 deg. curve lies above the 55 deg. curve. 


The line labelled IlLimit for Regular Reflection ll serves the same purpose as 
the similar line of Figs. 5 and 7.* In this case, for a given incident over- 
pressure, if the slope angle is smaller than the slope angle values on the line, 
Mach reflection will generally take place. 


3. 
Channeling. 
The last subject dealt with in the small charge studies was 


IIchannelingll in which a shock wave proceeds along a valley with parallel or con- 
verging walls. 
When this happens, the shock waves moving along the valley side 
walls reinforce each other, increasing pressure near the axis of the valley. 
Flat bottom, IIV 
II bottom, and converging valleys were studied, and in Ref. 12, 
simple relationships between pressure ratio and the slope angles of the valley 
side walls were developed. 
These relationships are shown in Fig. 13. 
The slope 
angles of importance are identified in the inserts on the figure. 


D. 
Techniques for Making Predictions 


1. 
Shock Front Pressures. 
The curves and plots relating overpressures on a 
slope and in channeling situations to various terrain and blast parameters all 
deal with a single topographic feature: a single change of slope from the terrain 
over which a shock wave is propagating, or a single elongated valley for channel- 
ing. 
But real terrain generally differs from this model, with rising slopes 
succeeding falling slopes, and valley and hill characteristics constantly changing. 
For this reason, in Ref. 13, techniques were developed for making predictions of 
at least peak overpressure in real terrain, and many of these techniques were in- 
corporated in Ref. 14 now in general use. 


Two initial steps must be taken to make predictions with shock waves encountering 


* Because the demarcation between regular and r1ach ref1ection is not sharp, the 
limit lines in Figs. 5 and 12 are not identical. The Mach reflection limit 
line (Fig. 5) lies above the regular reflection limit line (Fig. 12) at values 
of incident pressure less than about 50 kPa, and below that line at higher 
values of incident pressure. 
Values of reflected pressure from Figs. 5 and 12 
differ somewhat for these conditions. 
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Fig. 12. 
Overpressure at the Shock Front on a Rising Slope (Reflected Over- 
pressure) vs Peak Incident Overpressure for Shock Wave Undergoing 
Regular Reflection. 


Note: 
Numbers in parenthesis identify extreme (limit) slope 
angles along the "Limit for Regular Reflection" line. Other numbers 
on the figure identify the slope angles to which each curve applies. 
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slopes. First, a so-called IIprimary reference plane 
ll (PRP), must be establish- 
ed; second the effective slope an~le e (relative to the PRP) at a point for 
which predictions are to be made (a target) must be determined. 
The PRP is 
actually any convenient ground surface area over which shock wave characteris- 
tics can be determined, between ground zero and a target. 
(The area where 
a ~·1ach stem from an above surface burst first forms, for example, could consti- 
tute a PRP.) 


The effective slope angle e at a target may be determined most directly, as 
shown in the sketch, by drawing a line on a plan view of an area of interest 
whi ch connects GZ with the target. 


Line drawn perpendicular 
to contours at Target 
Target 


The angle between this line and a perpendicular to the contours of the target 
is the angle~. 
The modified effective slope angle e 
l 
, derived from a sectional 
view along this line is related to the nominal slope angle 8s by: 


tan e 
l = (tan 8s) (cos ~) 
(lb) 


The angle es can be calculated from this relationship, and used in Eq. 1a to cal- 
culate the effective slope angle 8. 


For most topographic conditions that might be encountered, however, the angle 


8 
1 is very close to 8, and can be used in its place (thus avoiding all calcula- 
tions except the initial determination of e 
l from the sectional view). 
It is 
only when both the angles es and ~ are large that the more complex procedure 
might be used. 
Table 1 illustrates these points. 


Once 8 is determined, Figs. 5 and 6 (based on vJhitham theory) or Figs. 7 and 
8 (based on small charge work) can be used to determine shock front pressure on 
a slope, using as incident pressure the value of peak pressure at the target 
location determined as if there were only flat terrain between GZ and the target. 


In channeling situations the prediction graph (Fig. 13) requires knowledge 
of the slope angles of the valley sides relative to the PRP, to be determined as 
is shown in the insets to the figure. If the direction of propagation of the 
shock wave is not within about 20 deg 
to the axis of the valley, the topo- 
graphic situation should be treated as a series of slopes with a shock wave 
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Tab 1 e 1 


Comparisons of .Modified Effective Slope Angle 8' With 
Effective Slope Angle 8, for Various 
Values of Slope Angle 8s and Azimuth Angle ~ 


Azimuth Angle ~ IModified Effective 
(deg.) 
Slope Angle 8 1 


( deg. ) 
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Effective Slope 
Angle 8 
( deg. ) 


9.4 
7.6 
5.0 
18.8 
15.2 
9.9 
28.0 
22.5 
14.5 
37.2 
29.5 
18.8 


traversing them at a large azimuth angle, instead of a valley in which channeling 
occurs. 


The characteristics of a shock wave after it has traversed a number of topo- 
graphic features is actually determined by its successive responses to each 
feature between GZ and the target area. It was shown in Ref. 25, however, (and 
in the current version of that document, Ref. 7) that for dete~mining the pressure 
at the front of a shock wave at a target, a "l ocal slope" approach could be used. 
In 
this approach, the shock front pressure is determined solely from the effec- 
tive slope angle at the target relative to the PRP, that is, as if the PRP was 
the slope immediately preceding that at the target. The general rationale for 
this approach is that effect of positive and negative slopes of equal angle 
essentially compensate for one another, so that pressure at the front of a shock 
wave moving over terrain features would be determined by the sum of the angles 
(positive and negative), of the slopes over which it had passed. 


The Whitham theory provides some support for this approach. 
In its appli- 
cation, the theory implies use of the approach, with e -- again either positive 
or negative -- being added to the Whitham theory function w to determine pressure 
at the shock front on a slope. 


It was noted in Ref. 14 that techniques which appear to be more rigorous (in 
which an attempt is made to follow a shock wave through a series of slope changes) 
would themselves only yield approximate values of shock front pressure. Super- 
imposed on any topographic effect would be the normal decrease in shock wave 
overpressure with distance from GZ, and no studies of these two phenomena occur- 
ring simultaneously have been made. 
Thus, the local slope approach, in addition 
to being simpler to use, may well be as accurate as any more complex approach. 


2. 
Pressures Behind the Shock Front. 
How to determine pressures behind 
the shock front after a wave has passed over a succession of slopes is a subject 
which has not been previously addressed. 
That the problem is complex can be seen 
from Figs. 9 and 10. Clearly, rising slope and falling slope effects are not com- 
pensating in many instances. 
At tiL of just under 0.1 mslm, for example, on a 
falling slope of 30 deg. the pressure ratios of shocks with incident pressures 
ranging from 14 to 173 kPa are all about 0.63, while on a 30 deg. rising slope, 
the pressure ratios at the same value of tiL range from about 1.3 to 2.2. Further- 
more, as is illustrated in Fig. 11, pulse shapes themselves differ at different 
points on the slope. 


In the absence of better information, a quasi-local slope approach was 
adopted for this program in an attempt to evolve a workable prediction method 
for predicting shock wave characteristics behind the shock front. 
As with pres- 
sures at the shock front, the effective slope angle at a target relative to the 
PRP is employed to determine which curve of Figs. 9 and 10 would be used. 
The 
value of a L, (which when multiplied by tiL from Figs. 9 and 10 determines 
actual time) is the distance to the target measured from the first slope change 
encountered, rather than from the slope change nearest the target. 
In addition, 
for stations other than the first encountered, no attempt is made to predict 
shock wave details between the time of shock front passage, and a time associated 
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with tIL values from Figs. 9 and 10, at which very rapid pressure ratio changes 
are no longer taking place. Further discussion of this method appears in a 
later section of this report in which all prediction techniques are evaluated. 


III. 
MODEL CONSIDERATIONS 


A. 
Model Characteristics and Scaling 


Since a major purpose of the entire program was to evaluate the adequacy 
of prediction methods in situations that might actually be encountered in the 
field, it was decided early in the program to base the model on a "real world" 
area that is of some interest to tactical planners. 
The prototype area chosen 
by BRL is shown in Fig. 14. 
The full-scale weapon yield of interest was 125 kt; 
the height of burst was 300 m. 
In order to enhance the terrain effects on blast, 
(the better to evaluate methods for predicting these effects) it was decided 
that the vertical scale of the model based on Fig. 14 would be increased by a 
factor of five. 


Scaling of shock wave phenomena in general depends on the fact that gravity 
has minimal effect of these phenomena. 
Ignoring other secondary effects (e.g., 
viscosity, real gas effects) and assuming that ambient atmospheric conditions over 
a model and a prototype are the same, the minimal effect of gravity implies that 
strict geometric scaling vli11 apply, that is, a geometrically scaled model \'1i11 
produce all desired phenomena without relative distortion.* In addition, energy 
W will scale as the cube of the length scale factor (the ratio between character- 
istic lengths in prototype and model) or -- equivalently -- the length scale 
factor between prototype and model equals the cube root of the ratio between ex- 
plosion energy yield in the two. 
In the tests with high explosives, the yield 
ratio is the same as the ratio between charge \'1eights. 


In the type of test being considered -- i.e., a burst over a surface -- 


all lenqths, including the burst height h, and distance from ground zero, d scale 
as t~ 1/3. 
At equal scaled distances from ground zero, (distance divided by W 1/3), 
shock pressures will be identical, and all times t, including pressure pulse 
duration t+ will scale as W1/3. 
In equation form: 


p = fl (d/W1/3) 


t/W 1/ 3 = f2 (d/W1/3) 
where d = distance from ground zero 


(3a) 


* Other explosion phenomena in which gravity is important cannot be scaled in so 
simple a fashion. 
For example, the cavity formed in a fluid by the burst of a 
high explosive charge is gravity dependent. 
The energy contained in the cavity 
scales as the fourth pO\,/er of the length scale factor. 
(The product of volume 
and depth is invarient.) But the energy of the explosive scales as the cube of 
the length scale factor. 
Thus, there will be distortion; the size of the cavity 
will change relative to the size of the charge that creates the cavity in two 
experiments using different charge sizes, all else being equal. 
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Fig. 14. Full Scale Terrain Contours Provided by BRL. 


If the atmospheric pressure (or densi ty) al so changes between two experiments, 
R.G. Sachs showed in Ref. 17 that the relationships still hold, in general, if 
shock pressure p is replaced by shock strength p/Po, and the measure of explo- 
sive yield ~~ is replaced by Po/l~ (nO\'/ more commonly W/Pp) where Po is the 
ambient pressure. 
Thus the pressure and time relationships become 


(p/Po) = f1 [d/(W/PO)1!3] 


and 
[t/(W/P )1/3] = f2 [d/(W/P )1/3J 
o 
0 


The quantity (W!po)1/3 has come to be known as the effective charge weight, 
since it serves the same function as actual charge weight in Eq. 3a. 


( 3b) 


Scaling from one nuclear burst condition to another, or from one chemical 
high explosive condition to another~ is straightforward with either Eq. 3a or 
3b. 
Scaling from nuclear explosions to HE explosions, while not as direct, has 
been accomplished by observing that about one half* of the energy in a nuclear 
explosion goes into forming blast waves. 
Thus, the weight of the high explo- 
sive charge equivalent to a nuclear burst with an energy yield of one kt 
(2 x 106 lb) would be about 106 lb or 4.54 x 105 kg. 


In the specific case being considered, a nuclear burst with a yield of 
125 kt is to be modelled by a high explosive charge weighing one lb or 454 gm. 
The length scale factor between full-scale and model-scale is the cube root of 
the ratio of the equivalent HE charge weights, viz (125 x 0.454 x 106/0.454)1/3, 
or 500. 
Thus the full-scale burst height of 300 m would become 0.6 m in the model. 


Based on these scaling considerations, the model derived from ng. 14 is 
shown in Fig. 15. 
All lengths (vertical as well as horizontal) are in metres; 
the zero elevation on the model would be the 200 m elevation full-scale. 
The 
numbers on F'ig. 15 show elevations from this datum, with vertical spacing between 
contours being 0.4 m. 
Elevat'ions of high points and low points are also shown. 


B. 
F1 at-Terrai n Reference Informati on 


The most frequently used measure of the effects of topography on a blast 
wave at a point on the ground surface is the change in peak overpressure of the 
incident shock wave, (the wave that would exist at the point if the terrain be- 
tween the point and ground zero were flat). Changes also occur in pressure 
behind the peak -- see, for example, Fig. 11 -- and BRL wished to determine 
whether methods available for predicting these changes were adequate. Thus, in- 
formation on the characteristics of the entire overpressure pulse over flat 
terrain was required. 


The relationship that wi'll be used in this report between shock front over- 
pressure on flat terrain and distance from ground zero from the burst of a 454 gm 


* This proportion differs for different weapon types, and heights of burst, and 
for different high explosives used to model a nuclear blast. It is, however, 
approximately true for the burst conditions of interest in this study. 
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Fig. 15. 
Contours for Terrain Model (200 m full scale taken as datum). 


charge 0.6 m above the surface, is shown in Fig. 16. 
Provided by BRL, this re- 
lationship is compatible with the height-of-burst (HOB) information in Refs. 19 
and 20 derived from charges ranging in weight from 454 gm to 114 kg.* 


At the pressure levels of interest in this report, the overpressure pulse 
from an explosion (nuclear or HE) has been described as a "modified exponential" 
pulse. Starting at a maximum at the shock front, pressure falls off with time 
very rapidly at first, then more slowly (as an exponential pulse does). After 
a certain time, however, dependent largely on explosive yield and distance from 
the explosion, the overpressure becomes zero. 
(An exponential pulse \A/ould only 
approach zero, thus the designation of modified exponential.) 


The time after the arrival of a shock front at any point on flat terrain 
that the overpressu~e falls to zero, is termed the positive pressure phase dura- 
tion, t+. 
The relationship between this parameter and ground range that will be 
used in this report is plotted on Fig. 16. 
For ground ranges greater than two 
metres it is the theorectical value from Von Neumann (Ref. 21); for smaller 
ground ranges, it is the curve given in Ref. 22 for bursts of hemispherical 
charges on a surface (which compares very well with that from Ref. 21 for ground 
ranges greater than two metres). Reported experimental values of durations show 
wider variation. It is a difficult parameter to measure because pulse overpressure 
approaches a zero value relatively slowly, and the zero crossing point is often 
difficult to identify. 


If both the overpressure and time after arrival of the pressure pulse at 
various distances from an explosion are normalized (by dividing pressures within 
the pulse by shock front overpressure, and time after shock front arrival by 
positive duration), it will be found, in general, that the higher overpressure 
pulses (those measured closest to the explosive source) are steeper at early 
times than lower overpressure pulses. 


A relationship between overpressure and normalized time that fairly well 
represents the pulse from either HE or nuclear explosions and for free-air, 
on-surface, or above-surface bursts is 


p(t) = p(o) (1 - 
T) e- UT 
(4) 


* These HOB curves contain information from a number of different chemical high 
explosives, the main ones being TNT and Pentolite. 
t~hile the characteristics 
of the explosives are by no means identical (TNT, for example, is oxygen de- 
ficient, and Pentolite releases some 10% more explosive energy than TNT), 
differences among the explosives appear to be masked by other experimental 
variables (including variations from batch to batch of the same explosive). 
At the scaled burst height used in this project, for example, mean ranges at 
which 30 to 50 kPa have been measured differ by more than 10% and individual 
measurements show a still wider variation. Therefore, the term high explosives 
as used in this report is meant to refer to the general class of high explo- 
sives with yields per unit weight about those of ,TNT or Pentolite. 
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where: 
p(t) = overpressure at any time t after passage of the shock front, 
thus 
p(o) = shock front overpressure, i.e., pressure at t = O. 
T = t/t+ 
t = time after arrival of the shock front 
t+ = positive phase duration 
a = a dimensionless variable dependent on peo) 


While this simple relationship is only an approximation of the pulse from 
an explosive source,* given normal variations from charge to charge, and the 
difficulty of determining a parameter as important as the duration t+, Eq. 3 
is deemed to be suitable for purposes of this report. 


The parameter a in Eq. 4 varies with shock front overpressure p(o) -- or 
incident overpressure Pi -- and the type of explosive source, principally 
whether HE or nuclear. Tabulated values in Ref. 23 for TNT (apparently derived 
from Ref. 24) indicate that: 


a '" 1 


0.04 p.O.55 
, 


for Pi < 140 kPa 


for Pi > 140 kPa 
(5) 


With a nuclear source, va'lues of a derived from Ref. 18 (as an approximation 
to the more complex relationship actually used in that report) are: 


0.14 p.O.51 
, 
for Pi < 70 kPa 


a '" 0.31 Pio. 32 
for Pi > 70 kPa 
(6) 


For values of Pi between about 40 and 700 kPa, values of a from Eq. 6 are 
greater than those from Eq. 5, which indicates that the nuclear pulse is gen- 
erally steeper immediately behind the front than the HE pulse. 
Below Pi '" 40kPa 


the HE pulse is initially steeper. 


Because of the minimum overpressure of interest (Pi'" 7 kPa) a gauge line 
extending 14 m from GZ was requireq for use in the model. 
Initially 15 
stati ons at fi xed di stances from G Z were to be used both on fl at terrai n and on 
topographic features. 
(T\I/o stations were added later, at BRL's request, to meet 
specific requirements.) 
To minimize the need for gauge location changes, the 
distances between adjacent stations of the standardized line were made symmetric 
about the center station (Sta. 8) -- that is, they are the same from either end 
of the line. 


~n Ref. 18, for example, the time history of a pulse from a nuclear weapon is 
given as: 


pet) = p(o) [ae-aT+ be-ST+ ce-YTJ [1 - TJ 


where a, b, c, a, B, 
and yare functions of incident overpressure. 


40 


Table 2 identifies the stations, and gives distances from ground zero as 
well as between adjacent stations. It also lists predictions of both peak 
flat terrain overpressures and positive overpressure pulse durations from Fig. 16 
for each station. 


In Fig. 17 are plots of flat-terrain pulse shapes for each station in 
Table 2. 
For clarity, time scales have been expanded at the higher overpressure 
stations because of the short pulse durations. 
To provide a sense of how sharply 
pulse durations actually do change, however, the last parts of pulses from sta- 
tions nearer ground zero have been superimposed on the plots for Sta. 4 and 
Sta. 6. 
It can be seen that the pulse at Sta. 4 is 7.5 times longer than that 
at Sta. 2, and the pulse at Sta. 6 twice as long as that at Sta. 4 (or 15 times 
longer than that at Sta. 2). 


C. 
Test Parameters - 
Burst and Gauge Locati ons 


The criteria used to establish burst and gauge placements reflected the 
basic purposes of the program: 
to determine whether currently available pre- 
diction techniques are adequate to deal with practical problems that might be 
encountered in the field. That simple statement led to three corollary criteria: 


1. 
Changes in blast parameters due to terrain should be significant. 
A change of 20% was adopted as the lower limit of significance. 


2. 
The terrain situations considered should include most (if not all) 
of the different types of situations that could be encountered. 
As a minimum, the situations should include effects on blast waves 
of both rising and falling slopes, successive slopes, contours at 
an angle to the shock front, location of a measuring point on a 
slope, and channeling situations. 


3. 
The overpressures considered should be those with which military 
planners must be concerned. 
This last criterion tended to put a lower limit on the incident 
blast pressures of interest. Again, a lower limit of 7 kPa(fI,,1 psi) 
was adopted. 


Plotted on Fi gs. 18, 19, and 20 is a set of seven G Z 1 ocati ons and ei ght 
gauge lines that meet these criteria. Three of the lines (A-2, B-2, and Ca-1) 
lie along essentially flat terrain; the others traverse various terrain features, 
and provide the desired variety of blast-terrain interactions. 
The dots along 
the gauge lines represent the locations of individual gauge stations identified 
in Table 2.* Note that gauge lines B-1, C-l, and C-4 have burst locations at 
both ends (designated as Band Ba, C and Ca' and C and Cb respectively). 


* Coordinates (including elevations) of each GZ and of each active station on 
each gauge 1 i ne are gi ven in Appendi x B. 
The x, y coordi nates refer to the co- 
ordinate system of Figs. 18, 19, and 20. 
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Table 2 
Standard Gauge Line Station 


Distances, Flat Terraln Pressures, and Duration 


Distance· 
Flat Terrain 
Positive 
Station 
from GZ 
Between 
Overpressure 
Pressure Pulse 


Number 
(m) 
Stations 
(kPa) 
(psi) 
Duration 
(m) 
(ms) 


1 
0 
13,800 
2000 
0.8 


2 
0.8 
2,070 
300 
0.2 
0.4 


3 
1.2 
690 
100 
0.6 
0.8 


4 
2.0 
210 
30 
1.5 
1.3 
5 
3.3 
70 
10 
2.3 
1.7 
6 
5.0 
33 
5 
3.0 


1.2 
7 
6.2 
24 
3.5 
3.2 


0.8 


8 
7.0 
21 
3 
3.4 


0.8 


9 
7.8 
17 
2.5 
3.6 
1.2 


10 
9.0 
14 
2 
3.8 
1.7 
11 
10.l 
10 
1.5 
4.0 


1.3 
12 
12.0 
9.0 
1.3 
4.2 


0.8 


13 
12.8 
8.0 
1.2 
4.2 


0.4 


14 
13.2 
7.5 
1.1 
4.3 


0.8 


15 
14.0 
7.0 
1.0 
4.3 
Added Stations 
4a 
2.6 
120 
17.4 
1.9 
5a 
4.1 
46 
6.7 
2.7 
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Fig. 20. Plan View, Gauge Lines 0-1 and 0-2, and Ground Zero O. 


Table 3 indicates which gauge stations along each gauge line are to be 
active during a test using the GZ shown at the top of the table. These stations 
correspond to the dots along the lines in Figs. 18, 19, and 20. 
The column on 
the extreme right of Table 3 indicates how many measurements equivalent to flat- 
terrain measurements would be made at each station location. 


In Figs. 21 through 25,profiles along each gauge line are plotted and 
elevation angles (the modified effective slope angles 8 
1 
) at each active station 
location are given. 
These elevation angles are tabulated in Table 4. 


D. 
Gauge and Recording System Considerations 


The basic purpose of the test program is to document the characteristics of 
the overpressure vs time pulses at ground level in the model, at gauge stations 
shown in Figs. 18 through 20. 
As will be seen in the section dealing with predic- 
tions, it is anticipated that some overpressure pulses will be marked by pressure 
spikes of relatively short duration -- about 0.1 ms at a number of stations, as 
short as 0.04 ms at one particular station. Whether the anticipated short dura- 
tion transients will be recorded and reported with an acceptable accuracy, depends 
critically on four elements: first, the response characteristics of the gauges 
or transducers that will sense the blast wave pressure; second, how these trans- 
ducers are to be used in the model; third, the system used to record the trans- 
ducer generated signals; and fourth, the characteristics of the system used to 
analyze the records. 
These elements are discussed, in turn, in the following 
material. 


1. 
Gauge Response. The response characteristics of the transducers planned 
for use by BRL are well suited to the requirements of this program. 
The PCB 
Piezotronics, Inc. series 100 transducers, to be used for lower pressure measure- 
ments, employ a quartz sensor, and various models have natural frequencies of 
between 250 and 500 kHz. 
The Susquehanna Instruments Model ST-4 transducers to 
be used for high pressure measurements, employ tourmaline sensors, and have a 
natural frequency of 1.5 MHz. 
The specifications for the PCB transducers quote 
their rise times as 1 ~s for the 500 kHz models, and 2 ~s for the 250 and 500 kHz 
models. 
The Susquehanna specifications merely state that " ... high speed reflec- 
tion shocks can be accurately measured to one microsecond or better." (This 
appears to be conservative; the rise time should be better than about 0.4 ~s.) 


2. 
Gauge Mounting. 
Given that the transducers have adequate response 
characteristics and that their function in this test program is to sense ground 
level pressures, the only requirements that apparently would need be imposed on 
their placement in the model would be that they be mounted flush with adjacent 
terrain, and that precautions be taken to avoid generating spurious signals 
(e.g., from gauge accelerations, cable noise, and the like). Unfortunately, 
doing so in this program might actually defeat the basic purpose of the program 
in some instances. 


The problem is purely a geometric one. 
The exposed area of the PCB trans- 
ducers (i.e., the diaphragms that transmit pressure generated forces to the sensing 
elements) have a diameter of about 0.55 cm. 
Shock wave velocities over the pressure 
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3 
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5.0 
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6.2 
7 


7.0 
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7.8 
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9.0 
10 
10.7 
11 
12 
12 
12.8 
13 
13.2 
14 


14.0 
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2.6 
4a 
4.1 
5a 
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per GZ 
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1 
1 
1 
1 


I 
1 
1 
1 


1 
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B 
B a 
B-1 B-2 
B-1 


Table 3 
Active Gauges for Each 
Gauge Line and Ground Zero 
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Gauge Line 


Station 
Number 


1 
2 


3 
4 
5 
6 
7 
8 
9 
10 


11 
12 


13 
14 
15 
4a 
5a 


A-I 
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0 
31. 0 
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, 


-24.2 
-24.2 


-15.6 
, 


0 
0 


I 


: 


a 


31.0 


Table 4 
Elevation (Modified Effective Slope) Angles at 


Stations Along Terrain Gauge Line* 
(degrees) 


B-1 
C-1 
C-2 
C-3 


0.7 
1.4 
1.4 


1.4 
-6.4 


- 0.7 


1.4 
-6.4 
1.4 
19.3 
-6.4 
41.2 
-36.3 
14.8 
0.7 
6.3 
-21.3 
14.0 
4.0 
28.8 
10.5 
14.0 
14.6 
5.3 
-16.7 


0 


a 
-9.1 


C-4 
0-1 


2.3 
2.3 


4.0 


16.7 
4.8 
-17.2 


4.0 


13.5 
0 


6.8 


-10.2 
-10.2 


a 
16.7 
-17.2 
* These angles are measured along each gauge line. They are not slope angles (i.e., elevation 
angles Os drawn normal to a contour as in Fig. 1), nor are they Effective Slope Angles 8 as 
defined in Eq. 1a and Fig. 2. 
They are the angle 0' given in Eq. 1b, and are very nearly 
equal to Os as shown in Table 1. 


0-2 


0.7 


0.7 


8.0 


8.0 


9.3 


-14.6 


range of interest (where the PCB transducer woul d be used), range from about 340 
to 600 m/s. 
Thus, the time required for the shock front to cross the sensitive 
area of the transducers range from about 9 to 16 ~s. (In the specific case in 
which a transient duration of 40 ~s is anticipated, the crossing time would be 
about 14 ~s, i.e., about 1/3 the duration of the transient.) 


If enough were known about the transducers (viz, their sensitivities with 
only part of their diaphragms subjected to pressure) and the velocity of the 
shock wave itself, the characteristics of a transient could be derived or recon- 
structed from a record. 
But there would always be some uncertainty, at least 
because an element used in the reconstruction -- shock front velocity -- depends 
on a quantity to be derived -- shock front pressure. 


For this reason, it is recommended that where short-time transients are ex- 
pected, an attempt be made to record head-on pressure at the shock front, from 
which actual shock-front pressure can be derived by use of the Rankine Hugonoit 
relationships for the shock front. 
An above-surface mounting block into which 
one of the PCB transducers could be placed in a head-on orientation would be 
small enough that any effects of its presence would not be detected at other 
stations. 
Such effects are well dissipated within about 10 block heights, and 
the minimum distance between adjacent stations is 40 cm (between Sta. 2 and 3, 
and 13 and 14, see Table 2). 
With duplicate test conditions, a head-on measure- 
ment might be made during one test to document the shock front, and a side-on 
measurement made during a second test to document the remainder of the pulse. 


3. 
Recording System. 
The recording system used to record the signals gen- 
erated by the transducers, must, of course, have response characteristics com- 
patible with quantity to be recorded. A typical magnetic tape system when 
operated in a standard FM record mode, has a maximum band width of about 0 to 
20 kHz (one db down), though an extended mode increases this to about 0 to 80 kHz 
(e.g., Honeywell Model 7610). A direct recording mode is also available whith 
changes the maximum band width to about 0.4 to 700 kHz. 
Other recording systems 
(e.g., cathode ray oscilloscopes) can increase the high frequency response still 
further. 


4. 
Data Analysis System. 
If the data analysis system incorporates some 
form of sampling at a fixed time rate, care must be taken that the system does 
not, in effect, act as a relatively low band-pass filter, which could very well 
modify (or totally eliminate) important pulse characteristics. If the recording 
system incorporates an analog mode, it is strongly recommended that this be 
played back and re-recorded for visual inspection at a proper speed. 
Visual 
examination is mandatory in any analysis system with inadequate (or borderline) 
sampling rate. 


IV. 
PREDICTIONS 


A. 
General Considerations 


Predicted shapes of the overpressure vs time pulses at those stations along 
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the 11 gauge lines -- 52 stations in all -- where overpressure pulses should 
be modi fi ed to a s i gni fi cant degree by the topography of the model in Fi g. 15, 
are given in Appendix A (along with the flat-terrain pulse at each station for 
comparison). 
For reasons to be discussed in the next section, where appropri- 
ate, predictions have been made of shock front pressures using both the Whitham 
theory curves (Figs. 5 and 6) and the curves based on small charge studies 
(Figs. 7 and 8). 
The shapes of later parts of the pulses were derived from 
Figs. 9 and 10~ Critical elements of the predicted pulse (e.g., shock front 
pressures, durations) are tabulated in Appendix B. 


In each set of traces (i,e., for stations along a particular gauge line 
using a particular GZ) predictions for the first station to experience any 
topographic effects are differentiated from those for stations further from GZ. 
(This is because predictions of pressures behind the shock front at these later 
stations were made with untried methods.) 
At each first station, the combination 
of shock front prediction made with the Whitham theory curves (Figs. 5 and 6), 
and later time predictions made with Figs. 9 and 10, is drawn as a solid line. 
At each more distant station that combination of predictions is drawn as a dashed 
line to indicate increased uncertainty. 
At all stations, however, where shock 
front predictions were made with the small charge curves of Figs. 7 and 8, the 
initial part of that pulse is shown as a dotted line. 


After finding shock front overpressures, a pressure ratio behind the front 
was found from Figs. 9 and 10 using a low value of tiL (between 0.1 and 0.3). 
Actual time was then calcu'lated by multiplying the value by L, the distance from 
the initial slope changes. 
Pressure in the incident pulse at that time was found 
from the plot of incident pressure vs time (Fig. 17), and actual pressure at 
that time was calculated by multiplying incident pressure by the pressure ratio 
from Figs. 9 and 10. 
Curves were then drawn connecting shock front pressures 
with the pressure-time point calculated as described. 
Simple curves or straight 
lines were employed to connect these two points except with stations on falling 
slopes for whi ch shock front pressures were deri ved from the Whitham theory 
curve, Fig. 6. 
In those cases, an initially sharply rising curve was employed 
to indicate the possibility of a rapid, early-time rise in pressure from shock 
front values. 
The reason for doing this is discussed more fully in the next 
section of this report. 


B. 
Discussion of Results 


The general natures of the predicted pressure-time pulses were as expected: 
shock front pressures increased on rising slopes, and decreased on falling 
slopes, and the greater the slope angle the greater the change from incident 
pressure. 
The characteristics of some of the prediction changes, however, are 
of some interest. 


For example, some of the pulses at stations on rising slopes are of very 
short duration. 
Sta 6 along line B-1 with GZ "B", and Sta 10 along line C-2 


with "C" (reproduced in Fi9. 26) are particularly noteworthy in that regard. 
At both stations, the spike is 23-24 kPa high based on Whitham theory (about 
11-13 kPa high based on small charge results). 
At Sta la, the spike duration 


* Special methods were used for Sta 11 and 12 of Line C-2. 
See Table B-2. 
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is about 0.08 ms, at Sta. 6, it is about 0.04 ms. 


Measurements on negative slopes may also display sharp pressure changes 
if Hhitham theory for pressure of the shock front holds. 
At stati ons such as 
No.8 along line A-I with GZ IIA II , No. 12 along line C-2 with GZ IIC II , and No. 15 
along line 0-2 with GZ "0 
11 (all reproduced Fig. 27), shock-front pressure pre- 
dictions from Whitham theory are 1/3 to 1/2 of those based on measurements with 
small charges. 
As discussed more fully in the next section of this report, 
this would imply very rapidly rising pressure close to the shock front. 


Finally, note should be made of the effects of channeling on the pulses 
along lines C-l and 0-1. 
Fig. 12 indicated that this valley would produce a 
pressure ratio of 1.25, that is, a 25% increase in overpressure. 
In the case 
of Sta. 10 along line C-l with GZ IIC II this was superimposed on a 70% increase 
(Whitham theory) due to rising slope effects to produce an overall increase of 
over 110% in shock front overpressure. (The overall increase, using small 
charge based shock front predictions, is about 80%.) 
In the case of Sta 15 
along line 0-1 with GZ 110 11 , the effect of channeling was simply to increase 
the incident pulse by 25%. 
Predicted pulses at both stations are reproduced 
in Fig. 28. 


V. 
APPRAISAL OF PREDICTION METHODS 


AND 


SOURCES OF NEW INFORMATION 


Among the more important findings of this effort to investigate and apply 
techniques for predicting terrain effects on blast, were: 
that there are unre- 
solved conflicts between two approaches, both of which had, at one time or an- 
other, been recommended for use in the field; that there could be some question 
about the accuracy of some of the available prediction curves; and that, for 
certain terrain situations, no techniques for predicting terrain effects had 
been developed. 
Since one of BRLls major purposes in undertaking the overall 
program was to uncover such deficiencies (if they existed) this phase of the 
program could be considered a success. 


In this section of the report the causes of these findings as they apply 
to predictions of pressures both at the shock front and behind the front are 
discussed first, followed by a brief discussion of ways in which needed addi- 
tional information could be acquired. 


A. 
Shock Front Predictions 


In Section II, two different sets of curves are given for predicting the 
effects of both rising and falling slopes: 
one from an approximate theory due 
to Whitham (Figs. 5 and 6); and one based on measurements with small charges 
(Figs. 7 and 8). 
The Whitham theory curves are more extensive than those 
from small charge tests, though both sets generally cover the overpressure 
range of interest. For the most part, Whitham theory predictions for rising 
slopes are higher, and for falling slopes are lower than the small charge pre- 
dictions. 
As can be seen in Figs. 29 and 30, for certain conditions, 
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differences in predictions of shock front pressure can be substantial. 


As already noted, many years ago it was recommended that the Whitham theory 
curves replace the small charge curves, because of both their theoretical basis, 
and the excellent agreement of predictions with shock tube measurements shown in 
Fig. 3. 
Upon review, however, whether there is adequate justification for this 
action is open to question, especially over the pressure region of interest 
(taken to be that for collateral damage, and damage to exposed material and equip- 
ment other than tanks, i.e., between about 10 and 100 kPa. 


On the one hand, while the agreement between Whitham theory and shock tube 
experiment shown in Fig. 3 is very good, the comparisons only extend down to 
shock strengths of about one, that is, to incident overpressure of about 100 kPa 
-- at the upper end of the range of interest. 
On the other hand, the measure- 
ments made with small charges were concentrated at overpressures between about 
10 and 200 kPa, that is, covering the entire range of interest. Thus, in this 
range, there is experimental evidence of shock front overpressure levels which 
do not agree with Whitham theory predictions. 


It is possible, however, that differences between the two sets of curves 
could be due to limitations inherent in the small charge test program, and 
thus that the conflicts between the two sources are only apparent rather than 
real. 
As can be seen in Fig. 10, and perhaps more graphically in Fig. 26, 


b~ast pulses on rising slopes can be characterized by pressure spikes at the 
shock front of very short duration. 
While the basic recording system used in 
the slnall charge studies was adequate to record such spikes, some rounding of 
the recorded pulse was inevitable (because of the time required for the shock 
front to pass over the sensitive area of the gauges used), and some subjective 
interpretation of the shock wave records was required (because each gauge dis- 
played different overshoot characteristics). 
In addition, the recording system, 
which was based on photographic recording of a single sweep of a cathode ray 
oscilloscope (eRO), could have masked some very short-time phenomena. 
Thus, re7 


ported shock front pressures on rising slopes from the small charge studies 
might be lower than they should have been, and curves based on this information 
would lead to low predictions of pressure. 


In a similar fashion, with falling slopes, the point at which maximum shock 
front pressure is attained can be very difficult to locate, if -- immediately 
after shock passage -- pressures rise rapidly. 
The finite time required for 
shock front passage over the gauge, and the method used to record the trace 
would both tend to mask the point. Thus, again because of gauge characteristics 
and recording system limitations, reported shock front pressures from the small 
charge studies might be higher than they should have been, and curves based on 
that information would lead to high predictions of pressure. 


B. Pressure Behind the Shock Front 


The prediction curves for pressures behind the shock front on falling 
slopes (Fig. 10) suffer from many of the same uncertainties as the curves for 
overpressures at the shock front. 
They were derived from tests with small 
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charges recorded in such a, way that subtle changes between incident pressures 
(also from measurements) and pressures on a slope might be difficult to find. 
The prediction curves for pressures on rising slopes (Fig. 9) do not, in gen- 
eral, have the same shortc:omings; they were derived from White's shock tube 
tests, which employed optical interferometric techniques instead of pressure 
transducers for determining shock strengths. 
The values at tiL = 0 for these 
curves, however, (that is, the values at the shock front) were derived directly 
from the small charge tests, and near that value of tiL, the curves should have 
increased uncertainty. Thus, there is some question about the overall accuracy 
of the falling slope curves, and about values near the shock front of the rising 
slope curves. 


In addition, as pointed out in Section II, there exists no simple technique 
to account for the effects of successive slopes on pressures behind the shock 
front, as exists for pressures at the shock front, where the local slope concept 
appears to work quite well. 
The quasi-local slope approach first used in this 
report to make predictions, is essentially an attempt to solve the problem by 
avoiding it, on the basic ground that further refinement is not likely to be 
warrented. 


What this approach does, for stations on slopes other than the first en- 
countered, is to delimit a portion of the pressure-time pulse immediately after 
passage of the shock front, during which predictions are just not attempted. 
The time length of this excluded portion of the pulse is based on the distance 
to a station in question from the first slope change that a shock wave encounters. 
Actual time after shock front passage is calculated by multiplying particular tiL 
values used with Figs. 9 and 10 by a length L. 
Use of a value of L based on the 
first slope change encountered in essence defines the portion of the pulse over 
which the influence of the first slope change is manifest. 
Other slope changes 
located between the first encountered and a particular measuring station will, 
at any given value of tiL) influence the pulse at shorter actual times, but no 
attempt is made to identify these influences. 


This technique has not been used before. 
Whether it yields reasonable 
results should be determined by the test program because of the very wide variety 
of topographic conditions built into the program. 


c. Acquisition of Additional Information 


It is anticipated that many of the conflicts and inaccuracies just discuss- 
ed can be resolved by the experimental phase of the overall program. 
As noted 
earlier, however, that phase could have some of the same experimental problems 
as those that led to the conflicts and inaccuracies in the first place, that is, 
transducer geometry and recording system frequency response might limit the 
accuracy of measurement of some of the short-time transients, and shock front 
pressures. 
(Transients are predicted to be as short as 40 ~s; time required for 
the shock front to pass over the sensitive area of a gauge will be from 9 to 16 ~s.) 
Ways of overcoming some of these experimental difficulties are discussed earlier 
in this report. 
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It is possible, however, that the experimental program could not resolve 
all conflicts to the accuracy desired, and additional work to do so would be 
required. A program to acquire the necessary new information could benefit 
both from lessons learned during the earlier program, and from advances in 
analytical capabilities made since their completion. 


For example, the earlier programs showed that information from two-dimen- 
sional 
analyses and tests could be applied to the basic three-dimensional 
situation inherent in the terrain effects problem, i.e., shock wave expanding 
spherically in free air encountering and propagating over assymetric topography. 
Therefore, much relevant information could be derived from two-dimensional 
approaches, namely, carefully thought out computer analyses, and shock tube 
programs. 
The latter would extend and improve White's and Bleakney's work, 
now some 30 years old (Refs. 11 and 26). 
Such programs could well be followed 
with additional high-explosive tests in either two-dimensional, or three-dimen- 
sional geometries, possibly with larger charges than will be employed in this 
program, as was done some years ago (Ref. 27). 
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Appendix A 


PREDICTED OVERPRESSURE vs TIME PULSES 


Included in this appendix are predicted pulses for all stations which 
should experience terrain effects along the gauge lines originating from 
the eight ground zeros shown in Figs. 18, 19, and 20, and tabulated in Table 3. 
For each ground zero, a sketch (from Figs. 21 through 24) and a key identifies 
the origin of the predicted curves. 
For the first station on each slope~pulses 


starting with Whitham theory predictions at the shock front are shown as solid 
lines; for all other stations, because of increased uncertainties, they are 
shown as dasRed lines. 
At all stations, pulses starting with shock front 
predictions based on small charge studies are shown as dotted lines. 
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Overpressures along Line A-I from a Burst over GZ IIA". 
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Fig. A-lb. 
Overpressures along Line A-I from a Burst over GZ IIAII (cont'd). 


71 


70 r 


60 I- 
1* 


60r 
ffii 
~ 
I 
@ 


= Incident 


50 8- 
50 r': 
- - - = Whi tham Theory 
4°k 
30 


40 1-\ 


Sta 6 


•••••••••• = Small Charge 


ro 
20 
0... 
""'" 
10 
OJ 
~ 
~ 
00 
Vl 
Vl 
1.0 
2.0 
3.0 
4.0 
OJ 
~ 
0 
-.....t 
~ 
N 
OJ 
> 
<:) 
25 


20 Po.. 
20 
Sta 8 


10 I- 
. .,---, 
101~ 
Sta 9 


o 
~ 
o 
1.0 
2.0 
3.0 
00 
1.0 
2.0 
3.0 
4.0 


Time (ms) 


Fig. A-2a. 
Overpressures along Line B-1 from a Burst over GZ "B". 


25 


Sta 10 
10L 
.. " 
..... 
Sta 11 


to 
0- 
~ 


(1) 
s- 
O' 
. , . 
. , . 
. , J=-e: 
, , 
0 
~ 
Vl 
0 
1.0 
2.0 
3.0 
4.0 
0 
1.0 
2.0 
3.0 
4.0 
5.0 
Vl 
'-I 
(1) 
W 
s-o. 
s- 
O) 
> 
0 
10 /-•••• ' .. 
10 


5 
Sta 12 
5 
Sta 13 


0' 
I. 
Ie 
7"'" 
o 
1.0 
2.0 
3.0 
4.0 
o I 
, 
• 
". 
I 
, 
. .,.....,., 
• 
I 


o 
1.0 
2.0 
3.0 
4.0 
5.0 


Time (ms) 


Fig. A-2b. 
Overpressures along Line B-1 from a Burst over GZ "B" (cont'd). 
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Overpressures along Line C-1 from a Burst over GZ "Ca" (cont'd). 
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Overpressures along Line C-3 from a Burst over GZ "C". 
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Appendix B 


TEST PROGRAM INFORMATION 


Tab 1 e B-1 - Coordinates of Ground Zeros and Pressure Gauges 


Table B-2 - Overpressure Pulses: Critical Pressures and Times 
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88 


00 
<.0 


Gauge line 
~ 
A-I 


Coordinate 
x 
y 
z 


Sta. No. 
1 
28.5 
26.3 
1.0 


2 
3 
4 
29.6 
24.7 
1.0 


5 
30.2 
23.5 
1.1 
6 
31.0 
22.3 
1. 95 
7 
31.8 
21.0 
1.85 
8 
32.2 
20.3 
1.4 
9 
32.7 
19.5 
1.20 I 
10 
33.2 
18.8 
1.0 


11 
34.0 
17.5 
1.0 


12 


13 


14 


15 
35.9 
14.5 
1.0 I 
Added Station 


5a 
30.6 
23.1 
1.6 


Table B-1 


Coordi nates* of Ground Zeros 


and Pressure Gauges 


A-2 
x 
y 
z 


28.5 
26.3 
1.0 


24.0 
22.3 
0.78 


22.0 
20.3 
0.67 
20.8 
19.3 
0.60 


B-1 
x 
y 
z 


46.8 
27.5 
0.80 


45.5 
26.8 
0.80 


43.6 
26.2 
0.85 


42.2 
25.8 
1.0 
41.0 
25.2 
1.5 


40.3 
25.0 
1.60 I 


39.6 
24.8 
1.25 


I 


38.5 
24.0 
0.80 I 
37.0 
23.9 
1.30 I 


35.5 
23.2 
1.5 
34.9 
22.9 
1.6 


33.5 
22.5 
1.7 


i 


* x,y coordinates are based on the coordinate system of Fig. 18,19, and 20. 
The z coordinate gives 
station elevation. 


<.0 
o 


Gauge line 
".. 
B-2 


Coordinate 
x 
y 
z 


Sta. No. 
1 
46.8 
27.5 
0.8 


2 
3 
4 
5 
45.5 
24.5 
0.84 


6 
7 
8 
44.2 
21. 2 
0.90 


9 


10 


11 


12 
13 
42.0 
16.2 
0.98 


14 


15 
41.5 
15.0 
1.0 


--------------- 
I 


Table B-1 (cont.) 


CoordinateS*of Ground Zeros 


and Pressure Gauges 


C1 


x 
y 
z 


33.2 
31.2 
0.88 


33.8 
\ 
30.3 
0.77 


35.3 
28.6 
0.54 


36.3 
I 27.2 
0.35 


37.5 
25.7 
0.55 


38.5 
24.0 
0.80 


40.8 
21.0 
0.90 


41.5 
20.0 
0.90 


C2 


x 
y 
z 


33.2 
31.2 
0.88 
I 


32.6 
29.7 
0.90 


I 
31.2 
24.6 
I 
0.98 


30.4 
22.7 
1.60 


29.8 
21.2 
2.20 


29.3 
20.2 
2.20 


28.5 
18.2 
1. 75 


* x, y coordinates are based on the coordinate system of Fig. 18, 19, and 20. 
The z coordinate gives 
station elevation. 


\.0 


Gauge line 
;0. C-3 
Coordinate 
x 
y 
z 


Sta. No. 
1 
33.2 
31.2 
0.88 


2 
32.5 
31.2 
0.90 


3 
4 
5 
30.0 
31.1 
0.95 


6 
7 
8 
26.5 
30.8 
1.27 
9 
10 
24.5 
30.7 
1.57 


11 
12 


13 
14 
15 
19.2 
30.4 
1. 78 


Table B-1 (cont.) 


Coordinates of Ground Zeros 


and Pressure Gauges 


C-4 


x 
y 
z 


33.2 
31.2 
0.88 


32.8 
32.2 
0.95 


30.4 
37.0 
1.32 


29.0 
39.0 
1. 90 


28.3 
40.4 
2.20 


27.8 
41.5 
2.20 


27.3 
42.3 
2.05 


26.7 
43.3 
1.80 


--- '----- 
------------ 


C -1 
a 
x 
y 
z 


41.5 
20.0 
0.90 


41.5 
19.2 
0.92 


41. 5 
18.0 
0.94 


41.5 
15.0 
1.0 


41.5 
12.2 
1.0 


41.5 
11.0 
1.0 


41.5 
8.0 
1.0 


41.5 
6.0 
1.0 


-------- 
~ 


* x, y coordinates are based on the coordinate system of Fig. 18, 19, and 20. 
The z coordinate gives 
station elevation. 


~ 
N 


Gauge line 
Coordinate 
Sta. No. 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


11 
12 


13 
14 
15 


,... 
0-1 


x 
y 


39.4 
37.5 
39.4 
36.9 


39.3 
35.5 
39.3 
34.4 
39.2 
32.8 


39.0 
28.6 


38.:5 
24.0 


Table 8-1 (cant.) 


Coordinate~ of Ground Zeros 


and Pressure Gauges 


0-2 


z 
x 
y 


0.95 
39.4 
37.5 


0.97 


1.10 
1.45 
1.10 
34.9 
40.0 


32.9 
40.8 


0.80 
31.2 
41.5 


0.80 
26.6 
43.9 


---- 


Added Stations 
4a 
39.3 
34.9 
1.3 
5a 
39.3 
33.6 
1.3 


z 
l 


I 
0.95 
! 
I 


1.03 


1.30 


1.68 


I 


1.62 


* x, y coordinates are based on the coordinate system of Fig. 18, 19, and 20. 
The z coordinate 
gives station elevation. 


Sta No. 
Pi 
t+ 


(kPa) 
(ms) 
1 
13,800 


2 
2,070 
0.20 


3 
690- 
0.60 


4 
210 
1.5 
4a 
120 
1.9 


5 
70 
2.3 
5a 
46 
2.7 


6 


\.0 
7 


33 
3.0 


24 
3.2 


21 
3.4 
17 
3.6 


14 
3.8 


10 
4.0 


w 
8 


Pi 


150 
150 
91 


120 
110 
47 
45 
42 
28 
5.2 
13 
10 


4.2 
11.5 
8.1 
5.2 
11.5 
6.0 


Pi 


Pi 


Table B-2 


Ovc:~pressu.re Pul ~~_~: _ f!"i!ic:~L ~ressures and Times 
----- 
-- --- 
-- 
"-- ---- 
---- -- -- -------- 


Line B -1 
GZ 'C' Line C-1 
-----, --. -- _._---_. 
, ... __ ._---- - -"-" ---,- 
pit) 
t 
p(o) 
p(t) 
t 


WT 
SC 
(kPa) 
(ms) 
(kPa) 
(kPa) 
(kPa) 
lms) 


Pi 
- 
- 


Pi 
I - 
- 


I 


- 
Pi 
I 


- 
- 
0.12 
Pi 


0.33 
0.40 
63 
53 
41 
0.04 
Pi 
- 
- 


0.96 
60(R) 
25 
0.40 


1.2 
1.0 
7.0 
12 
0.48 
48 
39 
23 
0.12 


1.4 
1.1 
10 
9.2 
0.64 
! 
28 
22 
10 
0.66 
30(C) 
25(C) 
14(C) 0.64 


23 
16 
6.6 
0.92 


---- 
--- 
--------~ --------- 
~----~ 
C-2 


GZ 'C' Line C-2 
-------- 


p(o) 
p(t) 
t 


iIT 
SC 
(kPa) 
(kPa) 
(kPa) 
(ms) 


Pi 
- 
- 


Pi 
- 
- 


, 
I 


I 


Pi 


42 
29 
18 
0.08 
30* 
20* 
8.0* 0.50 


9 


10 


11 


12 


13 


14 


15 


9.0 
4.2 
8.3 
4.2 
7.6 
4.3 
6.9 
4.3 


13 
11 
5.5 
1.1 
8.7* 
6.5* 
0.76 I 


Pi 


~: 


11 
10 
4.7 
1.2 


--~----~ 


Pi = Peak (shock fron~ pressure on flat terra1n 
p(a)= Zero-time (shock front) pressure on topographic 
feature 
t 
= Time after shock front passage 
p(t)= Pressure in pulse at time t 


12 
10 
4.1 
1.4 


Pi 
- 


t 
Positive overpressure phase duratl0n 
WT 
Pressures derived from Whitham Theory 


SC 
Pressures derived from small charge studies 
(C)= Channeling 


* These stations are at the foot of a very steep. conical feature with a very narrow base 
(diameter ~lm). Pressures at Station 11 drop rapidly to near incident values. 
Pressures at 
Station 12 should be close to incident values throughout the pulse. 


I 1 


- I 


',.0 ..,. 


Sta No. 


1 
2 


3 
4 
4a 


5 
5a 
6 


7 
8 
9 
10 


11 
12 


13 
14 


15 


"----"- 


Pi 
t+ 


(kPa) 
(ms) 


13,800 
2,070 
0.20 
690 
0.60 
210 
1.5 


120 
1.9 


70 
2.3 


46 
2.7 


33 
3.0 


24 
3.2 
21 
3.4 


17 
3.6 


14 
3.8 


10 
4.0 
9.0 
4.2 
8.3 
4.2 
7.6 
4.3 
6.9 
4.3 


Table B-2 (cont.) 


Ove!Q!'e5sure Pulses: Critical Pressure~ and Times 
-------------~-- 


GZ 'C' Line C-3 
GZ 'C' line C-4 
GZ '0' l ine 0-1 
GZ '0' line 0-2 I 
---- 
----------- 


p(o) 
p(t) 
t 
p(o) 
p(t) 


WT 
SC 
WT 
SC 
(kPa) 
(kPa) 
(kPa) 
(ms) 
(kPa) 
(kPa) 
(kPa) 


t 


ms) 


p(o) I 


WT 
SC 


~~) __ (kPa L 
I 


p(t) 
t 
p(o) 
p(t) 
t 
I 


WT 
SC 
(ms) I 
~kPa) 
(ms) 
(kPa) 
(kPa) 
(kPa) 


Pi 
Pi 
Pi 
Pi 
:: 
I 


Pi 


280 
250 
263 
0.03 
170 
160 
140 
0.09 


Pi 
80 
78 
60 
0.16 


23 
32 
29 
0.24 


15 
23 
19 
0.33 
42 
40 
35 
0.04 


Pi 
26 
24 
15 
0.56 
29 
25 
20 
0.24 


22 
18 
8.8 
0.96 
23 
19 
14 
0.20 I 
Pi 
22 
18 
12 
0.44 


12 
12 
8.0 
0.37 


2.2 
6.6 
3.8 
1.5 
1.9 
6.1 
2.9 
1.7 


2.0 
5.5 
3.7 
0.98 
8.7(C) 
I 
1.1 
5.0 
2.9 1.4 


-~-------- 
------"-.-.--~--------- 
._--- --- ------_ .. -- 
~-------------.- 
~: 
------"-:T------ 
P; = Peak (shock fron~ pressure on flat terrain 
t 
= Positive overpressure phase duration 
p(ol= Zero-time (shock front) pressure on topographic 
WT = Pressures derived from Whitham Theory 
feature 
SC = Pressures derived from small charge studies 
t 
= Time after shock front passage 
(C)= Channeling 


p(t)= Pressure in pulse at time t 


<..0 
(J1 


Table B-2 (cont.) 


Overpressure Pulses: 
Critical Pressures and Time* 


Sta No. 
t+ 


GZ 'Ba' Line B-1 


p{t) 


WT 
SC 
_________ LkPa) 
(ms) l_lkEal_l~~~ __ ~tk!'~l_{~s~ 


Pi 
p(o) 
t p;------ 


r- 
GZ C ' Line C-l 
GZ 'Cb' Line C-4 
p\o)-~ -p{tT--- t:- 
p(o) 
p(t) 
t 


WT 
SC 
WT 
SC 
(kPa) 
(kPa) 
(kPa) 
(ms) 
(kPa) 
(kPa) 
(kPa) 
(ms) 
6.9 
4.3 
2 
7.6 
4.3 
3 
8.3 
4.2 
12 
10 
5.5 
0.94 
1.8 
6.3 
4.0 
1.1 
4 
9.0 
4.2 
5 
10 
4.0 
Pi 
14 
12 
7.2 
0.73 


6 
14 
3.8 
5.0 
10 
6.1 
0.68 
18 
17 
11 
0.56 


7 
17 
3.6 
<1 
7.5 
7.8 
0.56 
6.4 
12 
11 
0.56 


8 
21 
3.4 
52(R) 
14 
0.96 
8.5 
16 
14 
0.47 


9 
24 
3.2 
62 
48 
18 
0.80 
10 
, 
33 
3.0 
23 
29 
19 
0.78 
34(C) 
38{C) 34{C) 
0.21 
9 
19 
19 
0.28 
111 
70 
2.3 
55 
60 
48 
0.33 
40 
52 
45 
0.11 


12 
210 
1.5 
Pi 
Pi 
13 
690 
0.6 
Pi 
Pi 
Pi 
14 
2,070 
0.2 


15 
13,800 
---______ 
1 ______ -'--- 


~: 
p. = Peak (shock fronij pressure on flat terrain 
p(a)= Zero-time (shock front) pressure on topographic 
feature 


t~ 
Positive overpressure phase duration 
WT 
Pressures derived from Whitham Theory 


* 


t 
= Time after shock front passage 
p(t)= Pressure in pulse at time t 


SC 
Pressures derived from small charge studies 


(C) = Channeling 


Note that for these ground zeros, values of Pi increase, and of t+ decrease as Station 
numbers increase. 


LIST OF SYMBOLS 


d 
Distance from ground zero 
K 
Variable in Whitham Theory 
L 
Distance from point of slope change 
M 
Mach number 
p 
Overpressure 
p. 
Peak (shock front) overpressure on flat terrain 
1 
p(t) 
Overpressure at time t 
p(o) 
Zero time (shock front) overpressure 


Po 
Ambient pressure 
t 
Time after passage of shock front 
t+ 
Positive phase overpressure duration 
w 
Whitham function 
W 
Charge weight or energy yield 


a 
A dimensionless variable in Eq. 4 


y 
Ratio of specific heats of a gas at constant pressure 
and volume 


8 
Effective slope angle 


8' 
Modified effective slope angle 


8s 
Slope angle; the angle of steepest assent 


~ 
Variable in Whitham Theory 


T 


Azimuth angle; the angle between the direction of 
shock wave propagation and a perpendicular to the 
slope contours 
t/t+ 
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Eq. 
3 
Eq. 
2 
Fig. 9,10 
Eq. 
2 
Eq. 
3 
Fig. 16 
Eq. 
4 
Eq. 
4 
Eq. 
3 
Fig. 9,10 
Fi g. 16 
Eq. 
2 
Eq. 
3 


Eq. 4,5,6 


Eq. 2 
Eq. la 
Eq. Ib 
Eq. la 
Eq. 
2 


Fig. 2 
Eq. 
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